Chaos in the random field Ising model 
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The sensitivity of the random field Ising model to small random perturbations of the quenched 
disorder is studied via exact ground states obtained with a maximum-flow algorithm. In one and 
two space dimensions we find a mild form of chaos, meaning that the overlap of the old, unper- 
turbed ground state and the new one is smaller than one, but extensive. In three dimensions the 
rearrangements are marginal (concentrated in the well defined domain walls) . Implications for finite 
temperature variations and experiments are discussed. 



The concept of chaos in disordered systems refers to 
the sensitivity of their equilibrium state (at finite tem- 
peratures) or ground state (at zero temperature) with 
respect to infinitesimal perturbations. In spin glasses 
H for instance it is well known that small changes of 
parameters like temperature or external field causes a 
complete rearrangement of the equilibrium configuration 
||,|]|. This has experimentally observable consequences 
like reinitialization of aging in temperature cycling ex- 
periments Q and has also been investigated in numerous 
theoretical works ||. 

A slight random variation of the quenched disorder has 
the very same effect on the ground state configurations. 
Although of similar origin, chaos with respect to temper- 
ature changes is harder to observe than chaos with re- 
spect to disorder changes ||, and the latter phenomena 
has been used to quantify spin glass chaos in numerical 
investigations 

This type of chaos has actually later been discovered 
in another, simpler random system, the directed poly- 
mer in a random medium [p|-pO|, which is equivalent to 
an dommain wall in a random bond ferromagnet. The 
interface displacement as a reaction to infinitesimal ran- 
dom changes of bond-strengths obeys particular scaling 
laws with exponents related to the well-known interface 
roughness exponent x • 

In this paper we consider the random field Ising model 
|jj and study, the first time to our knowledge, the sen- 
sitivity of its ground state with respect to small changes 
in the random field configurations. It turns out that 
the emerging picture is very reminiscent of chaos in spin 
glasses and random interfaces. This statement is quanti- 
fied by the following phenomenological picture [[|j9| . 

Consider a random Ising system defined, for instance, 
by the Hamiltonian 

(y) 1 

where Si — ±1 are Ising spins, (ij) indicates nearest 
neighbor pairs on a D-dimensional lattice of, say, linear 



size L, the Jij denote interaction strengths and hi local 
fields, both quenched random variables obeying some dis- 
tribution (continuous, in order to exclude ground state 
degeneracies). The case Jij Gaussian (with mean zero 
and variance one) and hi — is the spin glass (SG), 
the case Jij > and hi = is the random bond ferro- 
magnet (RBFM) and the case Jij — J and hi Gaussian 
(with mean zero and variance h r ) is the random field 
Ising model (RFIM). In order to study the sensitivity of 
the ground state of these systems with respect to small 
changes in the quenched disorder we can apply a random 
perturbation of amplitude 8 <C 1 to any of the quenched 
random variable. As a consequence the new ground state 
will differ from the old one. 

The RFIM ground state changes when the domain 
structure changes (for purely ferromagnetic states this 
argument does not work). One can estimate when the 
two ground states will be uncorrelated, beyond a length 
scale L* This can be found considering domain walls with 
an Imry-Ma |l|] type argument ||||: The energy Efi ip 
to flip droplets or domains or excitations of size L scales 
like L e , where 9 is the energy fluctuation exponent (9 
is denoted y in the SG context 

!)• The energy change 
due to the random perturbation E ran d scales like 5L d / 2 , 
where d = d s the fractal dimension of the droplet's sur- 
face in the SG case, d = D — 1 the interface dimension 
in the RBIM and d = D in the RFIM case. The decor- 
relation takes place when E ran d(L) > Efu p (L), i.e. for 
L > L* ~ S- 1 ^ with A = d/2 - y. In SG jargon L* is 
called the overlap length and A is denoted £, the chaos 
exponent [||J^]. 

Two remarks are in order: first, as pointed out in 
and jio) already for L < L* the ground state is slightly 
altered by the random perturbation. This is, however, an 
effect of the interplay between elastic energy and E ran d- 
This leads to displacements of the domain wall of size 
Ax ~ SL a with a = X + x, where y is the roughness 
exponent. The roughness exponent Jl^ ] and the energy 
fluctuation exponent 9 are related via 9 — 2x + D — 3 M . 

The second remark concerns the RFIM-case: for D < 
2 the concept of a macroscopic domain wall fails here, 
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FIG. 1. A ground state plus the perturbation-induced 

changes. The original spin-orientations are indicated in grey 

for ySj = +1 and white for Si — — 1, the flipped spins are 

indicated in black. L = 320, A = 2 and 5 — 0.1 (see text). 



and the above considerations can only be transfered cum 
grano salts. This means that they are sensible only for 
L <C £ ~ cxp(— h%/A), the typical size of domains in the 
2D RFIM " 
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In 3D the situation is different: The concept of a do- 
main wall is well defined and we get with the estimate for 
the roughness exponent x — 2/3 [[l2| and, consequently, 
9 = 4/3 the result A = 1/6, i.e. L* ~ (T 1 / 6 . The typical 
displacement of a domain wall thus is, as above, Ax ~ L a 
with a = 2/3 for L < L* and a = 5/6 for L > L*. 

If we take the above arguments serious for the RFIM 
in 2D the overlap length L* turns out to be formally in- 
finite (since with 9 = 1 one has A = 0), where one has 
of course to be careful due to logarithmic corrections to 
the energy Efu p . Thus the mechanism by which rear- 
rangements take place is due to the interplay between 
elastic energy and E ran d- Moreover, in the 2D RFIM, 
the typical displacement of domain walls should scale as 
Ax ~ SL for L <C £ since a = A+x = 1- As a conseqence 
the correlation or overlap between the old, unperturbed, 
ground state Si and the new one S-(<5) 



9 = ^E 5 ^) 



(2) 



behaves like 1 — q ~ L Q_1 and therefore q should be 
of order 0(1), depending on the probability with which 
domain wall displacements occur. 

In what follows we present results of exact ground state 
calculations for ID spin chains and for 2D systems. We 
use a random field distribution and a perturbation dis- 
tribution that have a constant probability density be- 
tween —A and A and —6/2 and 6/2, respectively and set 
Jij = 1 . Figure 1 shows an example of a large 2D ground 
state (L = 320) with the two spin orientations shown 
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FIG. 2. a) (Top) Scaling of the overlap parameter with 
random field strength for the ID spin chain. b) (Bot- 
tom) Scaling of the overlap in 2D for the system sizes L = 
40,80,160,240 and 320 and for 6 = 0.1. 



in white and grey, respectively, and the flipped spins in 
black. There are two features one should note. First, the 
size of the system is larger than the critical length scale 
needed for ground state breakup and the magnetization 
is practically zero. Second, the flipped spins form a num- 
ber of clusters of varying size, that seem to concentrate 
on the cluster boundaries of the original groundstate. 

Figure || shows what happens as one sweeps the RF- 
strength (A). In arbitrary dimensions, the limit A — > oo 
goes over to a site percolation problem, i.e. the local RF- 
orientation gives the spin state at a site. In that limit the 
overlap q is determined by the probability of the applied 
perturbation 6 to change the this orientation. For some- 
what smaller fields q gets smaller and in an apparently 
linear fashion as A changes. In the ID case the overlap 
is not sensitive to the system size above a certain thresh- 
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FIG. 3. Probability distributions of the overlap q for 
A = 1.8 for the system sizes L = 40 . . . 320. The in- 
set shows the standard deviations of the overlap pdf's for 
A = 1.4,1.6,1.8,2,2.4. 



old in A, below which the overlap quickly increases to 
unity again, which indicates a typical domain size. The 
overlap seems to become a (^-dependent constant in the 
thermodynamic limit and for A — ► 0. 

This ID behavior can be understood as follows. For 
simplicity let us assume that the first spin is fixed to be 
up, i.e. So = +1. Then the total random field energy 
at site n is given by H r = Xa=i hi m the unperturbed 
system and H' n = H n + A„ with A„ = Yl7=i m * ne 
perturbed one. If hi and Si are independently distributed 
variables with zero mean and variance h r — [hf] av and 
S r = [Sf] av , respectively, the variables H n and A„ are 
(for n ^> 1) Gaussian with mean zero and variance 
nh r and n8 r , respectively. The probability distribution 
P{H n , H' n ) is simply given by 

P{H ni H' n ) = J dA r P(H n )P{A n )S(H n + A n - H' n ) . 

Now the total RF fluctuations H n and H' n produce do- 
mains if their magnitude is large enough to overcome 
the ferromagnetic coupling: suppose that Si = +1 and 
Hi > — J for i = l,...n (i.e. a plus-domain), but 
H n+ i < — J; then S n +i will be flipped, i.e. S n +i = —1 
and a new (minus-) domain starts. For large enough typ- 
ical domain sizes the total RF fluctuations become large: 
one can neglect J and assume that only the sign of H n 
and H' n determines the ground state (note that this is 
different from the high field region h n 3> J, in which the 
local random fields hi dominate). Thus, the probability 
for S n and S' n being equal is given by 



p(S„ 



S' n ) = J dH n dH' n P{H n ,H' n ) ■ 9{H n ■ H' n ) , (3) 
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FIG. 4. Cluster size distributions of clusters of nipped 
spins for S = 0.1, L = 40, 160, 320, and A = 2. 



1 - ±8 r /hr + 0{5 2 ). For 
in which h 2 = A 2 /3 and 
we have 5 r /h r — 0.05 and 
i 0.97 agreeing roughly with 
-> in the limit L — > oo. 



where 9 is the step function. A straightforward calcu- 



lation yields p(S n — S' n ) = 
the data shown in Fig. la, 
5f. = 5 2 A 2 /\2 with 5 = 0.1, 
hence q = -l + 2p(S r = S' r ) p 
the numerical results for h r - 

The 2D behavior is depicted in Figure |for 5 = 0.1. 
The number of simulations is 10000 for L = 40 and 80, 
4000 for L = 160, 1000 for L = 240 and 500 for L = 320. 
The generic behavior of the overlap is as for the ID chain: 
q(A) is roughly linear until the regime of small fields 
(A < 2) after which it seems to saturate to a (5-dependent 
value q(S). The cross-overs (increase of q with decreasing 
A) are due to the ground state breakup mechanism. For 
small systems the ground state is ferromagnetic, except 
for a limited number of domains of the opposite spin ori- 
entation. The decrease in q is caused by the effect of the 
ground state becoming more and more uniform (mag- 
netization \m\ — » 1). Otherwise the behavior resembles 
strongly the ID case. 

The thermodynamic behavior of the overlap is also vis- 
ible in the statistics of overlap distributions. Figure |3| 
shows how the probability distribution P(q) of q behaves 
with varying system size and for A = 1.8 (the data is 
the same as presented in Fig. ||). For all systems P(q) 
is peaked at q — 1 but as L is increased a peak appears 
in the distribution, resembling a Gaussian. The inset of 
Fig. H shows the standard deviation 6q oiP(q) for varying 
A as a function of the system size L. Except for the by- 
now standard cross-over for small L and A we observe, 
that the width of the distribution decreases, which sig- 
nals that in the thermodynamic limit P(q) approaches a 
delta-function-like sharp one. The cross-over exponent c, 
defined with Sq ~ L~ c , seems to be exactly one (c = 1). 

The mechanism by which q is determined is illustrated 
in Figure ||. The size distribution of flipped clusters n(s) 
converges with L to a power-law, n ~ s -1 - 6 , with a cut- 
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FIG. 5. Dependence of the overlap q on the perturbation 
strength 5 for weak and strong magnitudes (A), L — 80. 



outcome is simply that of a random walk (Ax ~ L 1 / 2 ). 
In other words assuming that typical valleys in the en- 
ergy landscape are separated by an energy given by the 
energy fluctuation exponent gives completely different re- 
sults for temperature and ground state chaos than for 
random bond disorder. This discussion is intimately re- 
lated to coarsening and ageing in the RFIM; one should 
note that there are so far no simulation results that ad- 
dress these questions directly. 

ACKNOWLEDGMENTS 

This work has been performed within the Finnish- 
German cooperation project supported by the Academy 
of Finland and the DAAD. MA. would like to thank Eira 
Seppala for a version of the computer code used. 



off that depends very weakly if at all on L. This has 
to be so for the overlap not to diverge to zero in the 
thermodynamic limit, since one can write 1 — q as an 
integral over n(s): a L-dependent cut-off would imply 
that q would decrease continuously. 

Finally in Figure |^ we demonstrate that 1 — q ~ 5 for 
small 6. This follows from the scaling arguments pre- 
sented for 2D RFIM domain walls and the ID RF chain. 

In this paper we have considered the stability of the 
random field Ising model to small perturbations. Unlike 
in spin glasses, it turns out that the RFIM ground state 
shows a weak form of chaos, similar to directed polymers 
or random bond Ising model domain walls. The over- 
lap q attains its value from fluctuations of domain walls, 
in both ID and 2D. Thus the ground state stays almost 
intact. The ground state domains are robust against ex- 
ternal perturbations since, most likely, the field excess 
of a domain is extensive (^2 hi ~ V). For the RFIM in 
3D the prediction of the domain wall scaling argument 
is that q should converge to unity since the domain wall 
displacement exponent a here is 5/6: the displacement of 
a domain wall on large enough length scales is Ax ~ L a 
and therefore 1 — q oc L a ~ 1 — > 0. Moreover, in both limits 
h r /J—>0 and h r / J — > oc, i.e. deep in the ferromagnetic 
phase and deep in the paramagnetic phase, it is trivial 
that q — > 1. 

One would like to extend the argumentation to changes 
in temperature as is common for spin glasses and random 
bond -type directed polymers. In spin glasses chaos is in- 
timately linked to the non-equilibrium correlation length, 
which gives rise to measurable consequences in e.g. tem- 
perature cycling experiments that measure the out-of- 
phase susceptibility. Here, however, repeating the scaling 
argument of domain wall for temperature changes results 
in a displacement exponent which does not produce any 
extensive changes in the overlap. In 2D the predicted 
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